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METHODE D'ELABORATION DE CARTES DE PROBABILITE DE POSITION DE POINT D'IMPACT D'UN PUITS. 

&) Methode d'elaboration de cartes de probabilite de po- 
sition de point d'impact d'un puits. / 

Elle consiste a definir un point cible fixe et invariant, 
choisi a priori sur un emplacement initial E0 de la surface S, 
discretiser la surface S a I'aide d'un maillage compose de 
noeuds et de mailles, affecter au moins un vecteur d'incer- 
titude geometrique elementaire a chaque noeud du maillage 
de la surface S, determiner un nombre d'occurrences statis- 
tiquement significatif d' emplacements de la surface S en 
fonction des incertitudes geometriques qui I'affectent, proje- 
ter le point cible sur chaque occurrence d'emplacement de 
la surface pour en deduire un point d'impact, reporter Ten- 
semble des points d'impact sur Templacement E0 de la sur- 
face S en attribuant a chacun de ces points des 
coordonn^es surfaciques identiques a celles qu'il avait sur 
I'occurrence d'emplacement de la surface le contenant, de- 
finir sur la surface S une distribution statistique pour I'en- 
semble des points d'impact, determiner a partir de la 
distribution statistique une densite de probabilite" en tous les 
points de la surface S, ladite densite donnant la probabilite 
pour qu'un point quelconque de la surface S soit un point 
d'impact, et cartographier sur la surface S la density de pro- 
babilite* en courbes de niveaux. 
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Methode d'elaboration de cartes de probability de position de point 

<f impact d'un puits 

La pr£sente invention concerne une m&hode d'elaboration de 

5 cartes de probability de position de point d'impact d'un puits sur une surface 
incluse dans un milieu en trois dimensions et dont la position de chaque 
point est connue avec des incertitudes g^ometriques. 

Les incertitudes g^ometriques de position sont des vecteurs dont 
les composantes sont exprimSes en unites de longueur, par exemple en 

10 metres, dans les trois dimensions du milieu. 

Pour connaitre la structure g6ometrique du sous-sol, il est usuel, 
en exploration p&roliere notamment, de realiser ce qui est appete un bloc 
sismique 3D et l'ensemble correspondant de vitesses de sommation (stack). 
Les moyens d'obtention d'un bloc sismique 3D sont bien connus des 

15 sp^cialistes et ne seront pas decrits en detail. Ay ant obtenu un bloc sismique 
satisfaisant, on procede g6n6ralement & une etape de migration. L'opdration 
de migration a pour but de fournir h Tinterprdtateur une image du sous-sol 
qui se rapproche le plus possible d'une coupe g6ologique. 

Deux types de migration sont principalement utilises, & savoir la 

20 migration temps et la migration profondeur. La migration temps fournit une 
image acoustique. La migration profondeur fournit une image en 
profondeur du sous-sol aussi exacte que possible. Mais les deux types de 
migration n^cessitent au prealable, la definition d'un modele de vitesses, qui 
peut etre obtenu, par exemple comme proposd par LANDA (Geophysical 

25 Prospecting n° 86, pp 223-243, 1988) ou encore par WHITCOMBE (61* me 
SEG, Houston, 1991). 

En pratique, les modeles de vitesses sont entaches d'incertitudes, 
ce qui rend les r^sultats des migrations par la-meme impr£cis. Dans un 
article pr&ente a l'EAEG 1996 (Amsterdam 3-4 juin 1996), P. THORE et 

30 A. HAAS ont proposd une methode pratique et simple pour determiner les 
erreurs de migration dues a l'incertitude sur les vitesses, ainsi que des 
formules simplifies susceptibles d'etre utilises selon qu'on est en presence 
d'une migration profondeur ou d'une migration temps. Ces formules 
peuvent etre considerees comme une approximation de premier ordre pour 

35 la determination de la direction et du module d'un vecteur incertitude dans 
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le cas d f un champ de vitesses variable et pour autant que Ton soit dans le 
domaine de validity de la formule de DIX, c'est-£-dire avec des variations 
lat£rales uniformes de vitesses pour des reflecteurs plans et peu pentes. 

Dans le brevet FR-A-2 763 702, il est dScrit une methode 

5 d'elaboration de cartes de risques de positionnement d'un puits dans un 
milieu en fonction d'erreurs de migration dues k l'incertitude sur le champ 
de vitesses. La m&hode a pour objet d'effectuer autant de migrations d'un 
horizon sismique que Ton a de vitesses possibles et de noter les points 
d'impact d ! un puits de forage fixe sur chaque occurrence d'emplacement de 

io l'horizon migrd. 

Par ailleurs, la determination d'une trajectoire d'un puits de 
forage depend, d'une part, du choix d ! un point de depart et, d'autre part, du 
choix d'un point cible k atteindre dans le sous-sol. Classiquement, le point 
cible se situe sur une surface pr^sente dans le bloc sismique 3D, ladite 

15 surface constituant un horizon sismique. Etant donn6 que les emplacements 
des horizons sismiques sont affectes par des vecteurs d'incertitudes 
g6ometriques, le point cible ne peut avoir une position precise que sur une 
et une seule occurrence d'emplacement de l'horizon. En fait, chaque 
occurrence d'emplacement de l'horizon donne lieu & une nouvelle position 

20 du point cible et la reunion d'un grand nombre de ces positions offire une 
distribution des points d'impacts du puits de forage sur l'horizon. 

Dans la demande de brevet fran9ais numfro 99 02 088 du 19 
fevrier 1999, il est decrit une methode pour la determination d'une 
trajectoire optimale pour atteindre une cible floue dans un milieu a partir 

25 d'un point distant. La methode utilise les incertitudes geometriques affectant 
un volume pour determiner une trajectoire optimale d'un forage. Le volume 
est probabilisd a I'aide des incertitudes geometriques et la trajectoire 
optimale est celle qui traverse au mieux le volume probabilisg. 

La pr^sente invention a pour objet une m6thode permettant de 

30 dresser une carte de probability de position de points d'impact d'un puits de 
forage parfaitement d&ermind sur une surface incluse dans un milieu 
tridimensional d'axes x, y, z et dont Templacement est connu avec des 
incertitudes g6om6triques. En particulier, ladite surface est un horizon 
sismique inclus dans un bloc sismique deftni dans un rep£re tridimensionnel 

35 d'axes x, y, z et l'emplacement de l'horizon est entache d'incertitudes. 
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La methode selon l'invention concerne le risque de position- 
nement. Cependant, bien que son but soit de dresser des cartes de risques, 
elle ne realise pas pour cela des migrations successives selon diverses 
vitesses comrne dans FR-A-2 763 702. Plus pr£cis£ment, la methode 
5 remddie aux deficiences de la methode d^crite dans FR-A-2 763 702 en 
prenant en compte l'ensemble des causes d'incertitudes qui affectent 
Femplacement d'un horizon dans le sous-sol pour dresser une carte de 
probability de position de points d'impact. 

De plus, bien que la methode selon l'invention utilise les 
10 incertitudes g6om6triques, son but n'est pas de calculer une trajectoire 
optimale comme cela est decrit dans la demande de brevet fran^ais numSro 
99 02 088. 

La methode, selon l'invention, pour l'elaboration de cartes de 
probability de position de points d'impact d'un puits sur une surface S 
15 incluse dans un milieu en trois dimensions et dont I'emplacement est connu 
avec des incertitudes g6om6triques donndes, est caract6ris6e en ce qu'elle 
consiste k : 

- definir un point cible fixe et invariant, choisi a priori sur un 
emplacement initial E0 de la surface S, 

20 - discretiser la surface S a l'aide d'un maillage compose de nceuds 

et de mailles, 

- aflfecter au moins un vecteur d f incertitude geometrique 
eiementaire h chaque noeud du maillage de la surface S, 

- determiner un nombre d'occurrences statistiquement significatif 
25 d'emplacements de la surface S en fonction des incertitudes g6om6triques 

qui 1'affectent, 

- projeter le point cible sur chaque occurrence d f emplacement de 
la surface S pour en deduire un point d'impact, 

- reporter l'ensemble des points d'impact sur I'emplacement E0 
30 de la surface S en attribuant a chacun de ces points des coordonndes 

surfaciques identiques h celles qu'il avait sur l'occurrence d'emplacement de 
la surface le contenant, 

- definir sur la surface S une distribution statistique pour 
l'ensemble des points d'impact, 

35 - determiner h partir de la distribution statistique une densite de 
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probability en tous les points de la surface S, ladite density donnant la 
probability pour qu f un point quelconque de la surface S soit un point 
d'impact, et 

- cartographier sur l'emplacement EO de la surface S la densite de 
5 probability en courbes de niveaux. 

Selon une caracteristique de l'invention, les occurrences 
d'emplacements de la surface sont obtenues par tirages au sort des vecteurs 
d' incertitudes elementaires. 

Selon une autre caracteristique de l'invention, on attribue en 
10 chaque noeud du maillage de la surface plusieurs vecteurs d'incertitudes 
gyometriques elementaires et on calcule un vecteur incertitude global qui 
est la resultante desdits vecteurs d'incertitudes gyometriques elementaires et 
qui est affecte au noeud consider^. 

Selon une autre caracteristique, la surface S etant un horizon 
15 sismique defini par un pointe et une migration profondeur, on attribue & 
chaque noeud du maillage de l'horizon sismique trois vecteurs d'incertitudes 
geometriques elementaires, qui sont un vecteur d'incertitude sur le pointe de 
l'horizon, un vecteur d'incertitude sur la migration profondeur et un vecteur 
d'incertitude sur un calage sismique de l'horizon & au moins un forage du 
20 milieu, pour calculer le vecteur resultante desdits vecteurs. 

Selon une autre caracteristique de l'invention, en chaque noeud 
de l'horizon, les incertitudes gyom&riques 61ementaires varient en intensity 
sans varier en direction. 

Selon une autre caractyristique de l'invention, les variations de 
25 direction des incertitudes gyometriques elementaires en chaque noeud de 
l'horizon sont prydeterminyes. 

Selon une autre caracteristique de l'invention, en des noeuds 
voisins de la surface, les intensitys des vecteurs incertitudes rysultantes 
presentent des valeurs telles que lesdits noeuds ne se deplacent pas de fa9on 
30 independante les uns des autres lors de la determination de chaque 
occurrence d'emplacement de la surface. 

Selon une autre caracteristique de l'invention, lorsque la surface 
est traversye par au moins une cassure, les points de contact de la surface 
avec la cassure, ainsi que la partie du plan de cassure qui relie les morceaux 
35 de la surface entre eux sont solidaires de la surface lors de la determination 
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de chacune de ses occurrences d'emplacement. 

Selon une autre caracteristique de l'invention, la determination 
de la density de probability de presence d'un point d'impact quelconque sur 
la surface h partir de la distribution des points d'impact consiste a : 
5 - inscrire Tensemble des points d'impact dans un quadrilatere 

s'appuyant sur la surface S, 

- mailler le quadrilatere & l'aide de mailles ayant toutes la m€me 

aire, 

- attribuer a chaque maille une valeur de probability 6gale au 
10 rapport du nombre de points d'impact qui la traversent au nombre total de 

points d'impact inscrits dans le quadrilatere. 

Dans une mise en ceuvre avantageuse de la m&hode selon 
Tinvention : 

- on definit le point cible C comme un point de coordonnees Cx, 
15 Cy et Cz choisi a priori sur un emplacement initial d'une surface dans un 

milieu tridimensional, 

- on maille enticement la surface, 

- on attribue un ensemble (n-uplet) de vecteurs d'incertitudes 
g6ometriques 61ementaires h chaque nceud du maillage, 

20 - on proc£de k un tirage aleatoire des intensites des vecteurs 

incertitudes elementaires en chaque nceud de la surface puis on forme la 
r6sultante desdits vecteurs et on deplace chaque nceud selon le vecteur 
incertitude r£sultante qui lui est attribu£, chaque ensemble des nouvelles 
positions des noeuds definissant une nouvelle occurrence d'emplacement de 

25 la surface, 

- on projette verticalement le point cible sur chaque occurrence 
d'emplacement de la surface pour en deduire un nouveau point d'impact, 

- on transfere chaque point d'impact sur l'emplacement initial de 
la surface en lui attribuant, sur cet emplacement, des coordonnees 

30 surfaciques identiques k celles qu'il avait sur sa propre occurrence 
d'emplacement de surface, 

- on determine un quadrilatere circonscrit & l'ensemble des points 
d'impact et s'appuyant sur la surface, 

- on realise un maillage secondaire du quadrilatere a l'aide de 
35 mailles secondaries d'aires 6gales entre elles, 
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- on compte le nombre de points d'impact associds k chaque 
maille secondaire, operation qui produit une distribution sous la forme d'un 
histogramme dans les deux dimensions definissant la surface S, 

- on definit, k partir de Thistogramme, une densite de probability 
5 en deux dimensions par attribution a chaque maille secondare d'une 

probabilite 6gale au rapport du nombre de points d ! impact qui l'atteignent au 
nombre total de points d'impact contenus dans le quadrilatere, 

- on calcule la probability de presence d f un point d'impact en un 
lieu quelconque sur la surface en utilisant la densite de probability en deux 

10 dimensions, et 

- on trace des courbes d'isoprobabilite de points d'impact sur la 

surface. 

L'invention sera mieux comprise a la lecture de la description 
suivante d'un mode de ryalisation illustry par les dessins annexes sur 
15 lesquels : 

- la figure 1 est une reprysentation schymatique d'un bloc 
sismique 3D contenant une surface, figurant un horizon sismique, et un 
point cible, 

- la figure 2 represente un maillage d'un horizon sismique image 
20 du toit d'un ryservoir petrolier affecte par des vecteurs incertitudes en 

chaque nceud du maillage, 

- la figure 3 montre le toit et la base d'un reservoir petrolier dans 
sa totality, ledit ryservoir ytant sonde en quatre compartiments par trois 
failles, la figure 3a etant un agrandissement d'un contact du toit avec une 

25 faille, 

- la figure 4 montre le toit du ryservoir petrolier maille de la 
figure 2 traversy par deux failles et sur lequel on a indiquy un point cible et 
les points d'impact effectifs, 

- la figure 5 montre un quadrilatere maille circonscrit a plusieurs 
30 points d'impact et situe sur le toit du ryservoir petrolier, 

- la figure 6 montre des courbes d'isoprobabilite d'impact sur le 
toit du ryservoir pytrolier dyduites de la distribution des points d'impact de 
la figure 5. 

Sur la figure 1 est reprysente un bloc sismique tridimensionnel B 
35 d'axes x, y et z. Le bloc a yte obtenu par traitement conventionnel de 
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donnees de sismique reflexion acquises dans ce bloc et complete par une 
analyse effective par des interpretateurs. A l'interieur d'un bloc sismique, 
on distingue classiquement des horizons sismiques plus ou moins continus. 
Les horizons sont assimiles k des Iimites de couches de sediments et les 
discontinuity a des cassures ou des failles. 

A l'interieur du bloc sismique B, on a represente un seul horizon 
sismique, k savoir, la surface S figurant, par exemple, le toit d'un reservoir 
p&rolier. Cette surface est g6n6ralement affectee par des failles dont une 
seule F est representee sur la figure 1. 

La figure 1 montre egalement un point cible C choisi et fixe une 
fois pour toutes k l'interieur du bloc B. Le point cible C, consider^ comme 
invariant, est place sur un emplacement EO choisi a priori de la surface S. 
En r£alite, le point cible se trouve dans un voisinage immediat de la surface 
S puisque l'emplacement de cette derntere est connu de fa9on imprecise. 

La surface et les failles qui l'affectent sont discretisees par un 
maillage compose de noeuds et de mailles. Le maillage apparait comme une 
triangulation de la surface et il est d'autant plus fin que Tinteret petrolier est 
grand. En chaque noeud du maillage, la surface S se voit attribuer un vecteur 
d'incertitude geometrique global I qui varie en intensity et en direction 
(fig. 2). 

La position de chaque noeud de la surface est ainsi connue avec 
une incertitude geometrique globale I, laquelle est definie par un vecteur 
ayant trois composantes Ix, Iy et Iz selon les trois axes du repere du bloc 
sismique. 

Le milieu represente par le bloc B etant une portion du sous-sol 
et la surface S etant un horizon sismique defini par un points et une 
migration profondeur de donnees sismiques relatives a ladite portion, on 
attribue a chaque noeud du maillage de 1'horizon trois vecteurs d'incertitudes 
g^ometriques elementaires ou triplet, qui sont : une incertitude Ip sur le 
pointe des horizons, une incertitude Im sur la migration profondeur et une 
incertitude If sur le calage sismique par rapport a au moins un forage du 
sous-sol. Les incertitudes g6ometriques elementaires sont donn£es sous la 
forme de trois vecteurs, chaque vecteur etant accompagne de sa propre 
incertitude absolue. Par calage sismique on entend la mise en corres- 
pondance des horizons sismiques avec leurs emplacements observes dans 
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un puits de forage. 

Par exemple, le vecteur incertitude £16mentaire Ip sur le pointe 
possede une incertitude absolue dip. On consid£re alors dans la methode 
selon Invention, qu'une fonction de densit6 de probability est assoctee au 
5 vecteur Ip, laquelle est avantageusement une fonction de density gaussienne 
de moyenne nulle et d'ecart-type dip. 

En chaque noeud, on determine un module pour chaque vecteur 
d'incertitude elementaire par application d'une ra&hode de type Monte- 
Carlo qui associe un vecteur d'incertitude el6mentaire & un nombre tir£ au 

io sort compris entre 0 et 1 . 

En chaque noeud de la surface, on determine ainsi par tirage au 
sort les composantes de chaque vecteur raembre du triplet (Ip, Im, If) puis 
on forme la r6sultante I de ces vecteurs. On deplace ensuite chaque noeud 
dans le bloc B dans la direction du vecteur I et d'une quantity correspondant 

15 k son module. 

Quand on a effects les tirages au sort pour tous les noeuds 
compris dans la surface S, on obtient une occurrence d'emplacement de 
ladite surface S. Le nombre d'occurrences d'emplacement de la surface est 
consider^ comme statistiquement significatif lorsqu'on a determine au 

20 moins trente & cinquante occurrences d'emplacement ; de preference, le 
nombre d'occurrences d'emplacement est de 1'ordre de plusieurs centaines, 
voire de plusieurs milliers, cela dependant g£n£ralement du nombre de 
noeuds compris dans la surface S. 

Afin de ne pas ginivcr des surfaces non naturelles (trop 

25 froiss^es) lors des deplacements des noeuds, on attribue une certaine rigidity 
a la surface. Cette rigiditd est simulee par action sur les tirages au sort. 
Ainsi, au cours des tirages au sort des vecteurs incertitudes, notamment 
avec la methode de Monte-Carlo, on prend des valeurs tiroes au sort 
sensiblement proches les unes des autres pour des points voisins sur la 

30 surface. Cela signifie qu'en des noeuds voisins de la surface S, les intensites 
des vecteurs incertitudes r&ultantes pr&entent des valeurs telles que lesdits 
noeuds ne se deplacent pas de fapon independante lors de la determination 
de chaque occurrence d'emplacement de la surface S. 

La figure 3 illustre Taccrochage des horizons aux cassures qui les 

35 recoupent, le cas ech^ant. La figure montre un reservoir p£trolier limits 
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sup&ieurement par son toit ou surface M et inftrieurement par sa base ou 
surface N. Le reservoir est scindS en quatre compartiments par trois failles 
Fl, F2 et F3. Un grossissement du contact entre la faille Fl et le toit M (fig. 
3a), indique comment les maillages des horizons et des failles se relient 

5 entre eux par partage de noeuds communs. La portion de la faille Fl, 
marquee F 11, est consider^ comme faisant partie int£grante du toit du 
reservoir M lors de chaque d^placement du toit par tirage au sort des 
vecteurs incertitudes elementaires qui l'affectent. Par exemple, le noeud R, 
situe au contact du toit M avec la faille Fl, est commun aux maillages du 

10 toit M et de la faille Fl. De la sorte, les points de contact de la surface M 
avec la faille ou la cassure Fl, ainsi que la partie Fl 1 du plan de faille qui 
relie les morceaux de la surface M entre eux, sont solidaires de la surface M 
lors de la determination de chacune de ses occurrences d'emplacement. 

Dans une etape suivante, le point cible invariant C ayant pour 

15 coordonn^es Cx, Cy et Cz est projete de preference verticalement sur 
chaque occurrence d'emplacement de la surface S pour determiner un point 
d'impact Ensuite, l'ensemble des points d'impact ainsi determines est 
reporte sur 1'emplacement EO de la surface S en attribuant & chaque point 
d'impact, sur cet emplacement EO, des coordonn£es surfaciques identiques & 

20 celles qu'il avait sur l'occurrence d'emplacement de la surface S assoctee k 
ce point d'impact. Les coordonn£es surfaciques d'un point sur la surface S 
s'entendent comme des coordonndes par rapport a un couple d'axes portes 
par la surface S. 

La figure 4 repr&ente une vue de dessus d ! une surface maillee 
25 sur la base d'un maillage triangulaire. La surface est couple par deux failles 
Fl et F2. Au centre de la figure, un petit triangle represente le point cible C 
pris sur l'emplacement EO choisi a priori de la surface S. De part et d'autre 
du petit triangle, on observe un nuage de points d'impact figures par des 
petits ronds. Chaque petit rond represente la position d ! un point d'impact sur 
30 une occurrence d'emplacement de la surface, lequel emplacement est bien 
entendu different de l'emplacement EO de la surface choisi a priori comme 
support du nuage des points d'impact. L'ensemble des points d'impact de 
l'emplacement EO constitue une distribution statistique puisque chaque 
point d'impact est associd & une occurrence d'emplacement de la surface S. 
35 Le nombre de points d'impact est done statistiquement significatif. 
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A partir de la distribution des points d f impact, on definit une 
densite de probability en tout point de la surface S, ladite density perrnettant 
d'obtenir la probability pour qu'un point quelconque de ladite surface S soit 
un point d'impact. 

5 L'ensemble des points d f impact de la figure 4 situd autour du 

point cible C est inscrit dans un quadrilatfcre Q, par exemple un rectangle 
comme repr^sente sur la fig. 5. Le quadrilatere Q est maille a l'aide d'un 
maillage secondaire r^gulier, de pryfifrence avec des mailles carries ayant 
toutes la meme aire. Puis, on attribue & chaque maille du quadrilatere Q une 
10 valeur de probability d'impact 6gale au rapport du nombre de points 
d'impact qu'elle contient au nombre total de points d'impact inscrits dans le 
quadrilatere Q. 

La figure 5 montre un zoom sur un nuage de points d'impact. Le 
nuage est ici inscrit dans le quadrilatere Q qui supporte un maillage 

15 secondaire k base de mailles carries d'aire egale entre elles. Le comptage 
des points d'impact est realist maille par maille, par exemple la maille Ml 
contient 6 points d'impact, la maille M2 contient 2 points d'impact, etc. Les 
mailles etant toutes d'aire egale et les diverses positions de la surface S 
etant toutes 6quiprobables, on assimile le rapport du nombre de points 

20 d'impact contenus dans une maille donn£e au nombre total de points 
d'impact a la probability d'atteindre la maille donnye lors du forage. 
L'opyration de comptage du nombre de points d'impact associys h chaque 
maille du maillage secondaire du quadrilatere Q produit une distribution 
sous la forme d'un histogramme dans les deux dimensions qui definissent la 

25 surface S. A partir de l'histogramme ainsi produit, on determine une density 
de probability en deux dimensions. 

Dans une derntere etape, on cartographie sur 1'emplacement E0 
de la surface S la densite de probability en courbes d'isovaleurs de 
probability ou courbes d'isoprobability. 

30 La figure 6 reprysente le rysultat final a Tissue de l'application de 

la mythode. On observe le mSine nuage de points que celui de la figure 5 
ceint par des courbes d'isoprobability d'impact. La prycision de ces courbes 
est d'autant plus grande que le nombre de determination de points d'impact 
est elevy et que le maillage secondaire est serry. 
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REVENDICATIONS 



1. Mythode d'elaboration de cartes de probability de position de 

points d'impact d'un puits sur une surface S incluse dans un milieu en trois 
5 dimensions et dont 1'emplacement est connu avec des incertitudes 
g6ometriques, caract6ris6e en ce qu'elle consiste k : 

- definir un point cible fixe et invariant, choisi a priori sur un 
emplacement initial EO de la surface S, 

- discretiser la surface S k Taide d'un maillage compost de nceuds 
10 et de mailles, 

- affecter au moins un vecteur d'incertitude g^ometrique 
elementaire k chaque nceud du maillage de la surface S, 

- determiner un nombre d'occurrences statistiquement significatif 
d'emplacements de la surface S en fonction des incertitudes g£om6triques 

15 qui l'affectent, 

- projeter le point cible sur chaque occurrence d'emplacement de 
la surface pour en deduire un point d'impact, 

- reporter l'ensemble des points d'impact sur Templacement EO 
de la surface S en attribuant a chacun de ces points des coordonn£es 

20 surfaciques identiques k celles qu'il avait sur l'occurrence d'emplacement de 
la surface le contenant, 

- definir sur la surface S une distribution statistique pour 
Tensemble des points d'impact, 

- determiner k partir de la distribution statistique une densite de 
25 probability en tous les points de la surface S, ladite densite dormant la 

probability pour qu'un point quelconque de la surface S soit un point 
d'impact, et 

- cartographier sur la surface S la densite de probability en 
courbes de niveaux. 

30 2. Mythode selon la revendication 1, caracterisye en ce que les 

occurrences d'emplacements de la surface sont obtenues par tirages au sort 
des vecteurs d'incertitudes yiymentaires. 

3. Mythode selon la revendication 1, caracterisye en ce qu'on 

attribue en chaque noeud du maillage de la surface plusieurs vecteurs 
35 d'incertitudes gyomytriques elementaires et qu'on calcule un vecteur 
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incertitude global qui est la r£sultante desdits vecteurs d'incertitudes 
g£ometriques Elementaires et qui est affecte au noeud considErE. 

4. MEthode selon la revendication 3, caracterisEe en ce que la 
surface S etant un horizon sismique defini par un points et une migration 

5 profondeur, on attribue k chaque noeud du maillage de l'horizon sismique 
trois vecteurs d'incertitudes g£om6triques elementaires, qui sont un vecteur 
d'incertitude sur le pointe de Thorizon, un vecteur d'incertitude sur la 
migration profondeur et un vecteur d'incertitude sur un calage sismique de 
Thorizon k au moins un forage du milieu, pour calculer le vecteur rEsultante 
10 desdits vecteurs. 

5. Methode selon la revendication 4, caracterisEe en ce qu'en 
chaque nceud de Thorizon, les incertitudes gEomEtriques elementaires 
varient en intensite sans varier en direction. 

6. Methode selon la revendication 5, caractErisEe en ce que les 
is variations de direction des incertitudes gEomEtriques Elementaires en 

chaque noeud de Thorizon sont predetermines. 

7. Methode selon la revendication 4, caractErisEe en ce qu'en des 
noeuds voisins de la surface, les intensites des vecteurs incertitudes 
resultantes prEsentent des valeurs telles que lesdits noeuds ne se deplacent 

20 pas de fa9on independante les uns des autres lors de la determination de 
chaque occurrence d'emplacement de la surface. 

8. Methode selon la revendication 1, caractErisEe en ce que, lorsque 
la surface est traversEe par au moins une cassure, les points de contacts de 
la surface avec la cassure, ainsi que la partie du plan de cassure qui relie les 

25 morceaux de la surface entre eux sont solidaires de la surface lors de la 
determination de chacune de ses occurrences d'emplacement. 

9. Methode selon la revendication 1, caractErisEe en ce que la 
determination de la densitE de probability de presence d'un point d' impact 
quelconque sur la surface k partir de la distribution des points d' impact 

30 consiste a : 

- inscrire Tensemble des points d'impact dans un quadrilatere 
s'appuyant sur la surface S, 

- mailler le quadrilatere k Taide de mailles ayant toutes la meme 

aire, 

35 - attribuer k chaque maille une valeur de probability Egale au 
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rapport du nombre de points d'impact qui la traversent au nombre total de 
points d'impact inscrits dans le quadrilatSre. 
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